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1. Introduction 

'• In 1948, Casimir predicted the attraction between perfect metal parallel platedlldue 

. ^ to quantum fluctuations of the electromagnetic field. Recently, a multiscattering for- 

malism of the Casimir effect for the electromagnetic field has been presentecEEEIH 
The Casimir effect has some peculiarities. In particular, it is a non-pairwise 
interaction; the Casimir thermal force (the thermal part of the Casimir energy) 
between two isolating bodies is not necessarily monotonic in their separation, as 
seen in the sphere-plate and cylinder-plate cased^. In addition, for some geometries, 
intervals of negative entropy appear, as in the case of two parallel plates described 
by the Drude modeP or, as recently shown, in the interaction between a Drude 
model plate and spher^ and in the interaction between a perfect metal plate and 
sphere ■ . 

In this article we study the Casimir effect between two equal radii perfect metal 
spheres in the large and short distance limit. As a result, we find negative entropies 
in certain ranges of temperature and separation between the spheres. In addition, 
we find nonmonotonic behavior of the entropy with the separation while the force 
is attractive for all separations, making it appear as though the natural evolution 
of the system tends to increase the entropy in certain ranges of temperature and 
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separation. We discuss the thermodynamical meaning and consequences of negative 
entropies in Casimir effect. Similar results for spheres with more general dielectric 
models can be found in Ref. \TU\ and for Drude plates in Ref. [TT] . 

The remainder of the article is arranged as follows: In Sect. [2j we describe the 
multiscattering model used herein to obtain the Casimir energies and entropies for 
the two perfect metal spheres. In Sect. [3j we obtain the large and short separa- 
tion limit of the Casimir energy, In Sect. [4] we obtain the entropy of the system, 
analytically and numerically. In Sect. [5] we obtain the Casimir force. We discuss 
the thermodynamic consequences of these results in Sect. |6] Finally, we discuss the 
results obtained in the Conclusions. 



2. Multiscattering formalism of the Casimir energy 

To calculate the Casimir energy, entropy, and forces between two spheres, we em- 
ploy the multiscattering formalism for the electromagnetic fielcEEl The Casimir 
contribution to the Helmholtz free energy at any temperature T is given by 

oo ' 

E = kBTj2log\l-N{K„)\, (1) 

n=0 

where Kn — are the Matsubara frequencies and = 2JkBT thermal 
wavelength. The prime indicates that the zero Matsubara frequency contribution 
has height of 1/2. All the information regarding the system is described by the 
N matrix. For a system of two objects, this matrix is N = T1U12T2U21. is 
the T scattering matrix of the z— object, which accounts for all the geometrical 
information and electromagnetic properties of the object. Uij is the translation 
matrix of electromagnetic waves from object i to object j, which accounts for all 
information regarding the relative positions between the objects of the system. 

For a perfect metal sphere of radius R, the T matrix is diagonal in {ImP, I'm'P') 
spac^, with elements given by 



^EE __e . TT lIi+i{^R)-KRI^_.iKR) 

Expressions for the Vap matrices can be found in Refs. jSi and 1121 



3. Casimir energy 

3.1. Large separation limit 

To obtain the large separation limit of the Casimir energy, we need the dominant 
part of the T matrix in this limit. We define the adimensional frequency q by 
K ~ q/d, where d is the distance between the centre of the spheres. The main 



July 3, 2012 0:41 WSPC/INSTRUCTION FILE 

Casimir ■ Energy ' and ' Entropy 'between 'perfect ' metal ' spheres ' QFEXT 1 1 



Casimir Energy and Entropy between perfect metal Spheres 3 

contribution in the large separation limit comes from the lowest-order expansion 
term of the T matrix elements in 1 /d. The dominant contribution comes from the 
dipolar polarizabilities part of the T matrix, taking the form 

Using the property log \A\ — Trlog(A), and expanding Eq. ([!]) over ^, we obtain 

oo ' 

E = -ksTY^ Tr (N(Aj;i7i)) . (5) 

Using the translation matrices in a spherical vector multipole basi^ and the large 
separation approximation of the T matrix, the trace of the N matrix is obtained by 
straightforward calculus. Here we denote with a sub-index T the results valid for 
all temperatures, with a sub-index the results in the quantum limit (T 0), and 
with a sub-index cl the results in the classical limit (/i — > 0), which is equivalent 
to the high T limit. Carrying out the sum over Matsubara frequency, we obtain the 
Casimir contribution to the Helmholtz free energy for two spheres of equal radius 
R as 

^ _ HcR^ ze^^ ^ 
^ ~ ~2^2(e2^-l)''^ 

(2 (15 - 29z^ + 99z^) cosh(z) + 15 cosh(5z) + (-45 + 582^ + 182^) cosh(3z) 
-|-24z (6z2 - 5 + (5 + 3z^) cosh(2z)) sinh(z)) , (6) 

where z — d/Xr- We also define the adimensional Casimir energy as Ead{z) = 
^^Et- From this result, the quantum [T — > 0) and classical [h — > 0) limits with 
their first corrections are easily obtained as 

143;ici?^ ShcTT^R^ fkeTY 2288dncn'^R^ /ksT^^ 



^° 167rd7 27d \ he ) ^ 1575 [ he ) ' 

Ea = - ^^^J^ - ksT^ (15 + 30z + 29z2 + 18z3 + q^^) e'^^ , (g) 

Note that in Eq. Q, the first correction to zero temperature case is proportional 
to , contrary to the plate-sphere casd^Hl and to the cylinder-plate casd^, where a 
result proportional to were obtained for both systems. As corrections have the 
same negative sign of the Casimir energy in both limits, they describe an increase 
of the magnitude of the Casimir energy in high and low temperature limits. It is no 
longer the case if we study the energy for all temperatures. In fact, for some ranges 
of separation and temperature, thermal photons tend to reduce the Casimir energy 
between the spheres, as shown in the left plot of Fig. [TJ which is an indicator of the 
appearance of negative entropy in this system because of the negative slope of the 
energy curve. It is clear in Fig. [ijleft) that for some temperatures and distances, 
the Casimir energy between spheres is lower than in the zero temperature case. But 
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Fig. 1. The left and right figures arc the large {d,£ 2> R) and short (£ <g; R) distance limit of 
the Casimir energy between perfect metal spheres as a function of (left) and of (right) 
compared with the energy at zero temperature respectively. The dotted curve is the quantum limit, 
the dashed curve is the classical limit, and the solid curve is the asymptotic finite-temperature 
Casimir energy. Note that these curves are independent of the radius of the spheres. 



it is less evident the validity of Eq. ([7]). In fact, there is a tiny increase oi E/Eq at 
low until reaching a local maximun at ^ « 1.0388 of E/Eq « 1 + 10^**. It is 
not visible in Fig. [l] because of the difference of scales. 



3.2. Short distance limit, PFA 

To study the short distance limit of our problem, we apply the Proximity Force 
Approximation (PFA^JHEIEU We start with the general formula of Casimir energy 
between parallel plates for all temperature T as 



n=0 m=l •^'J 



(9) 



PFA is a valid approximation at short distances. Following Ref. UHl we define the 
distance between equal radius R spheres as h{9, (p) = £ + 2R{1 — cos{9)), where £ is 
the minimum distance between the surfaces of the spheres and p = Rsm{9), then 
we calculate the PFA energy as the short distance limit of the next integral 



E^J"^^ dip dpe|(£ + 2i?(l-cos(6l)), (10) 
Jo Jo 

where e\\ — -j^- Carrying out the integration over ip and applying the change 
t = 1 - cos(6'), 

= 27ri?2 f dt{t - l)ej{£ + 2Rt). (11) 

"'0 



t = and take the short distance limit {£ — > 0), leading to our final PFA result, 



To separate the divergence of the PFA energy, we apply the change of variable 
E^^^ = -2TTRt / duej{(\l + 2u)). (12) 
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The formula presented above is a valid approximation for all temperatures and 
material properties of the equal radius spheres at a short distance For perfect 
metal spheres we have Tr (i?i • R2) = 2, then the fc_L integral of Eq. (|9| can be 



carried out and the obtained result is applied to Eq. (12), obtaining 



= 4^ ^3 ■ (13) 

m— 1 n— 

The quanturrP and classical limits are 

^0 =-T777T^ ^ + ^ +"' (14) 



1440^2 72 \ he J 2 \ he ^ 

In Fig. 2l right), the energy in the short distance limit is shown as a function of 
In this hmit the energy grows monotonically with T, so we do not expect to find 
negative entropies in this limit. 



4. Casimir entropy 

In the canonical ensemble, the entropy is defined as S* = —drE. In the large 
separation limit, the Helmholtz free energy depends on the adimensional variable 
z = ~ 27r^^g^, so we can write the entropy as 

dz dE dkB dE 

and define the adimensional entropy as Sad{z) = j.'^^s S = —dzEadiz). From this 
result, the quantum (T 0) and classical {h 0) limits are easily obtained as 

^ ^ ^ lefcs^'i?' fkBTY ISSOAkBdTT^R^ ( keTV ^^^^ 
M \J^) 1575 \-h^) ' ^^^> 

= + fcs ^ (15 + 30z + 27 z^ + Uz^ + 9z' - ISz^) e'^^, (18) 

where z — So, the entropy is a growing function with temperature in both 
limits, but this is not the case for all temperatures, as we can observe in the left 
plot of Fig. [2] where a region of negative entropy and another region of negative 
slope of the entropy are observed. 

Because of the limit at low temperature of the entropy, we know that the entropy 
is positive for low T (not seen in Fig. [2] because it is small compared with S/Sd, 
but it can be observed in Fig. [S]), so there are three points where 5 = 0, including 
the origin. 

Negative entropy of the Casimir effect has already been obtained between Drude 
parallel plates in Ref. [7] and in Ref. 16, and more recently between a perfect metal 
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Fig. 2. The left and right figures are the large and short distance limit of the Casimir entropy 
between perfect metal spheres as a function of (left) and of (right) compared with the 
classical limit respectively. The dotted curve is the first temperature correction of the quantum 
limit entropy (the linear term of Eq. l |20[ l), the dashed curve is the classical limit, and the solid 
curve is the asymptotic finite-temperature Casimir entropy. Note the regions of negative entropy 
and negative slope of the entropy with the parameter ^ in the large distance limit dissappears 
in the short distance case. 



plate and sphere in Ref. W and between a Drude sphere and plate in Ref . [T7] and in 
Ref.m 

These results are only valid when the separation between the spheres is large 
compared with their radius, regardless of the radius of each one. 

In the Short distance limit, the entropy is directly obtained from Eq. (13) as 

S?'^ = E E (1 - 2-^^?^) ' (19) 

m— 1 n— 

and the high and low temperature hmits are 

° 36 \ /ic / 2 \ he J ' 

then the entropy is positive for all temperatures in the short distance limit. We 
can observe in Fig. [2j right) that the negative entropy interval dissappears at short 
distances and, at low temperature, the entropy is linear with T, contrary to the far 
distance limit, where a slope of was obtained. 



4.1. Numerical study at smaller separations 

As noted previously, asymptotic results are no longer valid when the separation be- 
tween the spheres becomes comparable to their radius. For this reason, a numerical 
study of entropy was performed for these cases. We computed Eq. ([!]) numerically 
for all temperatures from T = until reaching the classical limit for fixed ratio 
between the radius R and separation, i" = ^. 
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Fig. 3. The left figure is the log-log plot of the absolute value of the entropy divided by its classical 
limit for perfect metal spheres at constant r as a function of y~. Starting from the asymptotic 
solid gray curve (r — > 0), we increase r up to the dotted curve r = 0.4, observing a reduction of 
the interval of negative entropy. Dashed curve is the case r = 0.41, where the interval of negative 
entropy has disappeared. The solid black curve for r = 0.45 is shown and the PFA result for 
r = 0.45 is the dashed gray curve. The right figure is the log-log plot of the absolute value of the 
entropy divided by the PFA limit for perfect metal spheres at constant r as a function of too. 
In this case, the solid gray curve represents r = 0.05 instead the asymptotic far distance result. 



In Fig. [sjjleft) , the entropy of the system of two perfect metal spheres is plotted 
as a function of z = ^ for constant r — ^. The large and short separation limit 
result are shown too. We choose a log-log representation of the absolute value of the 
entropy divided by its corresponding classical limit. Therefore, zeros are observed 
as log divergences, and we can also observe the cases of negative entropy. Starting 
in the large separation regime, we observe an interval of negative entropy. As we 
increase r (reducing the separation between the spheres), the region of negative 
entropy tends to reduce until it disappears between r = 0.40 and r — 0.41. The 
power-law decay of the entropy at low temperatures (Eq. (17)) is observed as a 
linear decay of the curve at low z, and constant behavior in the high-temperature 
limit (Eq. (18)) is also reached in the computation. 

In the right plot of Fig. [3) we observe the ratio of the entropy with the PFA 
limit also as a function of z = ^ for constant r = ^. In this case we see that 
the numerical results tend to the PFA limit when spheres tend to contact on a 
non-trivial way, with a better convergence for high z than for low z. The power law 
decay of the entropy at low temperatures reduces from S oa for the large distance 
result (r — )■ 0), around S (x for the closest numerical studied case (r = 0.45) to 
S (X T for the short distance limit. This polynomial behavior S oc T"('') is observed 
in Fig. [3] as a linear decay of the curves at low z. As a consequence, the convergence 
of the entropy to the PFA result at low temperature will be reached at extremely 
short distances between the spheres. 



5. Casimir force 

In this section we calculate the Casimir force between spheres. In the large distance 
limit, the asymptotic Casimir force F = —ddE can be written in terms of the 
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Fig. 4. The left and right figures are the large and short distance limit of the Casimir force 
between perfect metal spheres as a function of (left) and of (right) compared with the 
energy at zero temperature respectively. The dotted curve is the quantum limit, the dashed curve 
is the classical limit, and the solid curve is the asymptotic finite-temperature Casimir energy. The 
nonmonotonic behavior of the force with temperature in the large distance limit results from the 
negative slope of the solid curve at constant separation and dissappoars in the short distance case. 



adimensional Casimir energy as 

Fadiz) ^'^^F^ lEadiz) - zd.Eadiz), (22) 

where z — In Fig.jiJ^left), the adimensional asymptotic force between the perfect 
metal spheres compared with the zero-temperature force is plotted as a function 
of for the large and short distance limits. Here, for constant temperature, we 
observe nonmonotonic behavior of the force with the adimensional parameter ^ 
in the large distance limit, which dissappears in the short distance case. 

It is easy to verify that the force behaves monotonically with separation, and 
the nonmonotonicity of the entropy with separation implies nonmonotonic behavior 
of the force with temperature (but not with the distance), because 

dT dTdd dd' ^ ' 

so the appearance of negative slopes of the entropy with separation implies non- 
monotonicity of the Casimir force with temperature, despite the attractive force for 
all separations and temperatures. 

As observed in Fig. |4j the force between the spheres is always attractive, but it is 
not always monotonic with temperature; asymptotically, for any given temperature, 
there exists a range of separations for which the force decreases with temperature. 
Then we have obtained a nonmonotonic behavior of the Casimir force with tem- 
perature between compact objects. Nonmonotonicities of the force between a plate 
and cylinder and between a plate and sphere were already obtained in Ref. [HI but in 
that case the nonmonotonicity already appears for the scalar field. Nonmonotonicity 
does not appear between spheres for the scalar field; this is a characteristic effect 
of the electromagnetic field, because cross-polarization terms of the Casimir energy 
are essential for this nonmonotonicity to appear. 
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In the Short distance limit, the force is obtained from the PFA energy given in 



Eq. ^ as 

F?'' - E E (1 + ■ (24) 



4£2 

m— 1 n— 



This result is shown in Fig. |4]J right), where er observe that the Casimir force is 
monotonous with the temperature in the short distance case. 



6. Thermodynamical Consequences 

In this section we discuss the thermodynamical consequences of the obtained results. 
In this article we have obtained the large and short separation limit of the Casimir 
energy and entropy for two perfect metal spheres. At any nonzero fixed temperature, 
we observe an interval of separations for which the entropy is negative, at zero 
temperature the entropy is always zero, while it is not the minimum of the entropy. 

In addition, the Casimir force is attractive for all separations and temperatures. 
So, we could naively think that we have possible violations of the second and third 
laws of thermodynamics, due to the existence of processes where the entropy of the 
system tends to increase and to the negative entropy intervals at finite temperature 
and distances, respectively. 

According to the Krein formulcli^li^ we know that the Helmholtz free energy of 
the electromagnetic field has three independent additive contributions: one is from 
the thermal bath, being proportional to the volume of the space^^. Another is the 
sum of contributions of objects immersed in the bath considered as isolated objects, 
each contribution being also a function of the volume and surface of each objectl^. 
The third is of geometrical nature, which we could call the Casimir part of the 
Helmholtz free energy and that we actually calculate in Eq. ([T]). 

Considering the whole system, the nonmonotonicity of the Casimir entropy with 
temperature is compensated by the contribution of the vacuum, because one scales 
with the global volume^ and the other with the separation between the spheres, so 
the behavior of the global entropy is dominated by the thermal bath contribution 
and the minimum is reached in the zero temperature limit, therefore there is not a 
violation of the third law. 

In addition, the second law of thermodynamics states that global entropy must 
increase for any process, but only in closed systems. As we are working in the 
canonical ensemble, we are implicitly assuming that there exists an external reser- 
voir which keeps our system at constant temperature, so the system is not isolated 
and its entropy can increase or decrease without violation of the second law. In the 
canonical ensemble, the condition equivalent to the second law is that the global 
Helmholtz free energy must decrease for any process, and this is true for the studied 
system. Therefore, the appearance of nonmonotonic entropy behavior in the canon- 
ical ensemble just implies nonmonotonic behavior of the force with temperature, as 
seen in Eq. (23). 
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7. Conclusions 

In this article we have obtained the large and short separation limit of the Casimir 
energy and entropy for two perfect metal spheres. At any nonzero fixed temper- 
ature, we observe an interval of separations for which entropy is negative, while 
at zero temperature the entropy is always zero. We showed numerically that there 
exists a minimum separation between the spheres for which the negative entropy 
interval disappear. In Sect. |6j having into account the complete thermodynamical 
system, that the system is described by the canonical ensemble, and with the help 
of the Krein formulcEffiH and Weyl formulcf^^, we have demonstrated that there are 
not violations of second and third laws respectively. Therefore, the appearance of 
nonmonotonic entropy behavior in the canonical ensemble just implies nonmono- 



tonic behavior of the Casimir force with temperature, as seen in Eq. (23 1, which 
is attractive for all separations and temperatures. While this proceeding was on 
review process, a related work was publishecP^l. 
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